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We present a measurement of D°-D° mixing parameters using the ratios of lifetimes extracted from a sample 
of D° mesons produced through the process D* + — > D°tt + , that decay to K~tv + , K~ K + , or 7r _ 7r + . The 
Cabibbo-suppressed modes K~ K + and tv~-k + are compared to the Cabibbo-favored mode K~ir + to obtain 
a measurement of yep, which in the limit of CP conservation corresponds to the mixing parameter y. The 
analysis is based on a data sample of 384 fb _1 collected by the BaBAR detector at the PEP-II asymmetric- 
energy e+e~ collider. We obtain y C p = [1.24 ± 0.39(stat) ± 0.13(syst)]%, which is evidence of D°-D° 
mixing at the 3cr level, and AY = [-0.26 ± 0.36(stat) ± 0.08(syst)]%, where AY constrains possible CP 
violation. Combining this result with a previous BABAR measurement of yep obtained from a separate sample 
of D° -> K~K+ events, we obtain y C p = [1.03 ± 0.33(stat) ± 0.19(syst)]%. 
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AY: 



Several recent studies have shown evidence for mixing in 
the D°-D° system at the 1% level fl, H H . This is consis- 
tent with Standard Model (SM) expectations |4S] and provides 
strong constraints for new physics models 15J. One conse- 
quence of D°-D° mixing is that the D° decay time distribu- 
tion can be different for decays to different CP eigenstates |@] . 
An observation of CP violation in D°-D° mixing with the 
present experimental sensitivity would provide evidence for 
physics beyond the SM Q]. We present a measurement of this 
lifetime difference and results of a search for evidence of CP 
violation in D°-D° mixing. 

The two neutral D mass eigenstates \Di) and \D 2 ) can be 
represented as 



ID,) = p\D°)+q\D ) 
\D 2 ) = p\D°)-q\D ) 



(1) 



where \p\ 2 + \q\ 2 = 1. We characterize the rate of D°-D° 
mixing with the parameters x = Am/T and y = AT/2T, 
where Am = mi — m 2 and AT = Ti — T 2 are the differ- 
ences between the mass and width eigenvalues of the states 
in Eq. (JTJ, respectively, and V = (Tj + r 2 )/2 is the average 
width. If either x or y is non-zero, mixing will occur. 

The effects of CP violation in D°-D° mixing can be pa- 
rameterized in terms of the quantities 



and 



iff = arg 



q_Af_ 

P A f 



(2) 



where A f = (f\H D \D°) (A f = (f\H D \D )) is the ampli- 
tude for D° (D ) to decay into a final state /, and H d is the 
Hamiltonian for the decay. A value of r m 7^ 1 would indicate 
CP violation in mixing. A non-zero value of iff would in- 
dicate CP violation in the interference between mixing and 
decay. Within the SM, CP violation in decay is expected 
to be small in the D°-D° system |8[] and is considered else- 
where n. 

D a -D° mixing will alter the decay time distribution of D° 
and D° mesons that decay into final states of specific CP. 
To a good approximation, these decay time distributions can 
be treated as exponential with effective lifetimes Tu h and rZ h , 



given 



by 1 



'hh 



'h h 



tk-k [1 +r m (y cos cpf - x sin iff)] 1 
tkh [l + r" 1 (y cos iff + x sin ip f )] 



(3) 



where tk-k is the lifetime for the Cabibbo-favored decays 
D° -> K--K+ and D° -> K+vr", and r,+ (r^) is the life- 
time for the Cabibbo-suppressed decays of the D° (D°) into 
CP-even final states (such as K~K + and These ef- 

fective lifetimes can be combined into the quantities yep and 
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[Thh 



(4) 



\Thh) 



r + 

'hh 



where (r hh ) = (r,+ + r hh )/2 and A T = (r,+ - r,^)/(r / 
T hh ). Both yep and AY are zero if there is no D°-D° mixing. 
In the limit of CP conservation, yep = y and AY = 0, with 
the convention that cos tpj > 0. 

We measure the D° lifetime in the three different D° de- 
cay modes, K~ir + , K~K + , and 7r~7r + . We use D° mesons 
coming from D* + — > D°tt + decays fioll ; the requirement of 
a D* + parent strongly suppresses the backgrounds. We use 
the charge of the £>* ± to split the K~K + and samples 
into those originating from D° and from D° mesons for mea- 
suring the CP-violating parameters. To avoid potential bias, 
we finalize our data selection criteria, the procedures for fit- 
ting and for extracting the statistical limits, and determine the 
systematic errors, prior to examining the mixing results. 

Most systematic errors related to signal events are expected 
to cancel in the lifetime ratios. Background events can contain 
effects that differ in each decay mode, making them difficult 
to characterize. Therefore, the event selection is chosen to 
produce very pure samples. The decay time distribution of 
signal candidates is fit to an exponential convolved with a res- 
olution function that uses event-by-event decay time errors. 
The decay-time resolution parameters are allowed to vary in 
the fit. Residual background components are modeled using 
Monte Carlo (MC) simulated events and control samples ob- 
tained from the data. 

We use 384 fb -1 of e + e~ colliding-beam data recorded 
near ^/s = 10.6 GeV with the BaBar detector HH at the PEP- 
II asymmetric-energy storage rings. We begin by reconstruct- 
ing candidate D° decays into the final states K~ir + , ir~ir + , 
and K~ K + . We require tracks to satisfy particle identifica- 
tion criteria based upon dE/dx ionization energy loss and 
Cherenkov angle measurements. We fit pairs of tracks with 
the appropriate mass hypotheses to a common vertex. We re- 
quire the invariant mass of a candidate track pair to be within 
the range 1.78-1.94 GeV/c 2 . To further reduce backgrounds, 
we require the helicity angle 8h, defined as the angle between 
the positively charged track in the D° rest frame and the D° 
direction in the lab frame, to satisfy | cos 6*^ | < 0.7. This is 
particularly helpful for rejecting combinatorial background, 
especially in the ir~7r + mode. 

We reconstruct D* + candidates by combining a D° candi- 
date with a slow pion track (denoted 7T+), requiring them to 
originate from a common vertex constrained to the e + e~ in- 
teraction region. We require the 7r+ momentum to be greater 
than 0.1 GeV/c in the laboratory frame and less than 0.45 
GeV/c in the e + e~ center-of-mass (CM) frame. We perform 
a vertex-constrained combined fit to the D° production and 
decay vertices, requiring the x 2 -based probability, P(x 2 ), to 
be at least 0.1%. The decay time t and its estimated uncer- 
tainty a t for each D° candidate are determined by this fit. 
We reject slow electrons that fake 7r+ candidates using dE/dx 
measurements in the tracking volume and further veto any 



5 



7r+ candidate that may have originated from a reconstructed 
gamma conversion or 7r° Dalitz decay. 

To reduce combinatorial backgrounds from D°'s produced 
via B-meson decay we require each D° to have a momen- 
tum in the CM frame greater than 2.5 GcV/c. We also require 
—2 < t < 4 ps and a t < 0.5 ps. The most probable value of 
at for signal events is 0.16 ps. For cases where multiple D*^ 
candidates in an event share one or more tracks, we retain only 
the candidate with the highest P(x 2 )- 

The distribution of the difference in the reconstructed D* + 
and D° masses (5m) peaks near 145.4 MeV/c 2 . Backgrounds 
are suppressed by discarding D* + candidates with a value 
of 5m deviating more than 0.8 McV/c 2 from the peak. In- 
variant mass distributions for the selected D° candidates are 
shown in Fig.Q] For the lifetime fit, we only use events within 
15MeV/c 2 of the D° signal peak (shaded regions in Fig 1); 
the event yields and purity within this signal region are also 
given. 
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FIG. 1: The reconstructed D° mass distributions for the three D' 
samples, within ±0.8 MeV/ c 2 of the peak of 5m. The shaded region 
indicates the events used in the lifetime fit. (The structures appear- 
ing above 1.92 GeV/c 2 in the K~ K + decay mode, and below 1.81 
GeV/c 2 inthe 7r ty decay mode, are mainly due to candidates with 
misidentified kaons or pions.) Also shown are the yield and purity 
of the three D° samples as calculated inside the ±15 MeV/c 2 mass 
window. 



The D° lifetime is determined from an unbinned maximum 
likelihood fit to the reconstructed decay time and its estimated 
error for events in the signal region. The fit is performed si- 
multaneously to all five decay samples (D° — * K~ K + ; D° 
-> K+ K~; D° -> 7T- 7T+; D° -> tt+ tt"; D° -> R--K+ 
and D Q — > K + ir~ combined). The D° candidates in the 
signal region can be divided into three components: D° sig- 
nal events, combinatorial background, and mis-reconstructed 
charm events. Each component is described by its own prob- 
ability density function (PDF) which also depends upon the 
D° or D° decay mode. 

The measured decay-time distribution of signal events is 
described by an exponential convolved with a resolution func- 
tion. The resolution function is the sum of three Gaussian 
functions with widths proportional to cr t . The three Gaussian 



functions share a common mean which is allowed to be offset 
from zero in order to take detector misalignment effects into 
account. The effect of the offset is studied as part of the cross- 
checks and taken into account as a systematic uncertainty. The 
resolution function parameters are all permitted to vary in the 
fit. Up to an overall scale factor in the width, the resolution 
function is observed to have the same shape for all modes, in- 
cluding the offset. To account for the small (1.5%) differences 
in width, we introduce two parameters Sk-k+ an d S n - n + 
to scale the overall width of the K~K + and tt~it + resolu- 
tion functions relative to the width of the K ~n + resolution 
function. All other resolution function parameters are shared 
among the different modes and are determined by a simulta- 
neous fit to all modes together. 

About 0.4% of the D° signal in the K~K + and 7r~7r+ 
modes consists of a correctly reconstructed D° combined with 
an unrelated 7r s ; this is estimated from MC and verified in 
data. These candidates have the same resolution and lifetime 
behavior as those from correctly reconstructed D* + decays, 
but about half of them will be tagged as the wrong flavor. We 
therefore include a 0.2% component in the signal PDF that 
uses the lifetime of the opposite flavor state. 

The decay-time distribution of the combinatorial back- 
ground is described by a sum of a Gaussian and a modified 
Gaussian with a power-law tail to account for a small num- 
ber of events with large reconstructed lifetimes. The means of 
these functions are allowed to float in the fit. Each of the three 
decay modes has its own shape for the combinatorial back- 
ground. These shapes are determined from fits to the events in 
the sideband region defined by 1.89 < M hh < 1.92GeV/c 2 
and 0.151 < 5m < 0.159 GcV/c 2 . We determine the amount 
of combinatorial background using MC samples scaled to the 
same luminosity as the data, modeling all known, relevant 
physics processes. The fraction of combinatorial background 
in the K'-n + mode is estimated to be (0.032 ± 0.003)%, in 
the K~K + mode (0.16 ± 0.02)%, and in the 7r _ 7r+ mode 
(1.8 ± 0.2)%. The uncertainties are determined by compar- 
ing data and MC events in the (Mhh, 5m) sideband where the 
combinatorial background is dominant. 

Mis-reconstructed charm background events have one or 
more of the charm decay products either not reconstructed or 
reconstructed with the wrong particle hypothesis. Most are 
D° mesons from a D* + — ► D°ir s decay with a correctly 
reconstructed tt s . For the K~tt + mode, most of the charm 
background is semileptonic decays D° — * K~ l + v with the 
charged lepton misidentified as a pion. The semileptonic de- 
cays also contribute to the K~K + final state, but the domi- 
nant contribution is from D° — > K~ir + ir° in which the 7r° 
is not reconstructed and the 7r + is misidentified as a kaon. 
There is also a small contribution from D + — > K~ir + ir + de- 
cays. In the 7r~7r + mode, the charm background is almost 
exclusively due to mis-reconstructed D° — > K~ir + decays in 
which the kaon has been misidentified as a pion. The decay- 
time distributions of the charm backgrounds are described by 
an exponential convolved with a Gaussian. The parameters 
are fixed to values obtained in a fit to MC events. The fraction 
of charm background events in the signal region is estimated 
from MC simulation and crosschecked by comparing data 
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and MC events in a (Mhh, Sm) sideband region defined by 
1.78 < M hh < 1.80 GeV/c 2 and 0.14 < 8m < 0.16 GeV/c 2 , 
where the charm background is the dominant contribution. 
We estimate the charm background to be (0.009 ± 0.002)% 
of events in the signal region for K~tt + , (0.2 ±0.1)% for 
K-K+, and (0.15 ± 0.15)% for n-n + . 

The results of the lifetime fits are shown in Fig. [2] The fit- 
ted D° lifetime t Ki is found to be 409.33 ± 0.70 (stat) fs, 
consistent with the world-average lifetime lfl2Tl . From the fit 
results we calculate yep and AY for the K~K + mode, the 
n~n + mode, and the two modes combined, taking into ac- 
count any correlations between the fitted lifetimes. The dom- 
inant correlation of 1 1 % arises primarily because the decay- 
time resolution offset is shared between the decay modes. The 
yep and AY results are listed in TableQ] The combined result 
is obtained by fitting the data with common lifetimes for the 
K~K + and tt~tt + modes, and assuming the same value of 
ipf for the K~K + and ir~n + decay modes. 



TABLE I: The mixing parameters extracted from the fit to data, 
where the first error is statistical and the second is systematic. 

Sample yep AY 

K~K + (1.60 ±0.46 ±0.17)% (-0.40 ± 0.44 ± 0.12)% 

7I--7T+ (0.46 ±0.65 ±0.25)% ( 0.05 ± 0.64 ± 0.32)% 

Combined (1.24 ± 0.39 ± 0.13)% (-0.26 ± 0.36 ± 0, 



Various cross-checks have been performed to ensure that 
the fit is unbiased and the assumptions in the fit model are 
well-founded. An offset in the resolution function is measured 
in the fit to be —4.75 ± 0.51 fs. This offset was seen in our re- 
cent K~tt + mixing analysis 01 and has also been observed 
in other BaBar measurements of charm decays. Because we 
measure ratios of lifetimes, the presence of a common offset 
has minimal impact on the values yep and AY. However, 
differences in the offset between the three decay modes, or 
between the D° and D°, could introduce a bias. No resolu- 
tion offset is found in the MC samples. However, we are able 
to introduce offsets in the fits to the MC sample of up to twice 
the size of the offset in data by misaligning the Silicon Vertex 
Tracker (SVT). In all cases the offsets are found to be consis- 
tent between all modes. 

The fitting procedure has been validated with generic MC 
samples weighted to the luminosity of the data sample and 
with dedicated signal MC samples. The signal efficiency is 
found to be independent of the true decay time and the fitted 
lifetimes are consistent with the generated value. 

The assumption that the resolution function is the same for 
all decay modes except for a scale factor is tested by fitting 
each sample independently. This gives mixing parameters and 
resolution offsets consistent with the nominal fit, but with sig- 
nificantly larger statistical uncertainties. The lifetime has also 
been extracted in independent fits to the flavor-separated sam- 
ples of D° -> K-TT+ and D° -> K+ir' decays. The fitted 
lifetimes and resolution functions in these two samples are 
consistent with each other. 
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FIG. 2: Decay time distribution in the data samples with the com- 
bined fit overlaid. The top left plot is the tagged K ~7r + sample, the 
middle plots are the D* + (left) and D*~ (right) tagged K~ K + sam- 
ples, and the bottom plots are the tagged -k~it + samples. The shaded 
and black distributions represent the charm and combinatorial back- 
ground in the fit, respectively. The normalized residuals for each fit 
are shown as a separate histogram for each sample. The top right plot 
shows a summary of the measured lifetimes. 



To cross-check the effect of the resolution offset, we per- 
formed further studies by dividing the data sample into sub- 
samples with different sensitivities to detector effects and fit- 
ting each subsample independently. Besides the D* tagged 
samples used for this mixing measurement, we also use a con- 
trol sample of D° — > K~tt + decays where the D° is not re- 
quired to come from a D* decay. This untagged sample has 
about five times as many D° decays as the D* tagged samples 
combined, allowing us to divide the sample more finely. The 
quantities used to divide the data into subsamples for these 
tests include the run period, the azimuthal and polar angle of 
the D° meson, and the orientation of the D° decay plane with 
respect to the X- Y (bending) plane of the detector. In all of the 
variables mentioned, the resolution offset is observed to have 
a large variation (typically between — lOfs and Ofs), but the 
fitted lifetimes are consistent among samples. Furthermore, 
the weighted average of the mixing parameters from the sub- 
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divided data samples is in almost all cases nearly identical to 
that obtained by fitting the full data sample with one common 
lifetime and resolution function as described previously. The 
largest variation is observed with the polar angle of the D° 
meson in the laboratory frame, where decays perpendicular 
to the beam line are found to have almost no resolution off- 
set, while decays into the forward region of the detector have 
a large offset. Since the acceptance for D° — > K~K + de- 
cays is lower in the forward region than for D° — > K~ir + or 



TABLE II: Summary of systematic uncertainties on yep and AY, 
separately for K~ K + and tt^tv + and averaged over the two CP 
modes, in percent. 



decays, the polar angle dependence in the offset 



could potentially introduce a different average offset for each 
of the three modes. This is accounted for in the systematic 
errors. 

The systematic uncertainties on the mixing parameters are 
small since most uncertainties in the lifetimes cancel in the 
ratios. We have considered variations in the signal and back- 
ground fit models, changes to the event selection and detector 
effects that could introduce biases in the lifetime. Table [TT] 
summarizes the various systematic uncertainties. The evalua- 
tion of each of these is described below. The systematic un- 
certainty on yep and AY averaged over the two CP modes is 
occasionally smaller than the individual uncertainties because 
of anti-correlations. 

We vary the signal PDF shape, and the size and position 
of the signal region. As part of the PDF shape variations, we 
perform a fit without a resolution offset. The effect of the 
polar angle dependence in the resolution offset is evaluated 
by performing the fit with separate, floating offsets in seven 
bins of polar angle, but sharing all other resolution parame- 
ters and lifetimes across all polar angle bins. The difference 
in the mixing parameters between this fit and the nominal fit 
is found to be small (< 0.02%). The largest systematic con- 
tribution to Dqp (0.12%) is due to widening the signal region 
mass interval from 15 to 25 MeV/c 2 . The choice of signal re- 
gion determines the level of mis-reconstructed signal events 
included in the fit. 

The mis-reconstructed charm background is a very small 
component in the lifetime fit and is determined using MC 
events. Varying the charm background fraction (depending 
on the mode) and the effective lifetime, both within their as- 
sociated uncertainties, yields a minor contribution to the sys- 
tematic uncertainty. 

Because of the high purity, the results have little sensitivity 
to the modeling of the combinatorial background, except in 
the 7r~7r + mode where varying the fraction of combinatorial 
background by 10% yields a systematic uncertainty in ygp of 
0.14%. We also alter the fit procedure by using a different 
sideband region and by substituting the MC decay time distri- 
bution for that obtained from fitting the data. Neither variation 
contributes a large systematic uncertainty. 

We have studied the effect of varying the event selection 
criteria, which could potentially affect the lifetime measure- 
ment. Changing the treatment of events where multiple D* + 
candidates share one or more tracks (either keeping all of them 
or throwing them all out) has little effect, while changing the 
upper bound on the decay time uncertainty from 0.5 to 0.4 ps 
yields the largest systematic uncertainty on y^ P of 0.172%. 
As with the D° mass window, the choice of the a t range af- 
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VCP (%) 




a AY (%) 




Systematic 


K'K+ 


TT TT + 


Av. 


K~K + 




Av. 


Signal model 


0.130 


0.059 


0.085 


0.072 


0.265 


0.062 


Charm bkg. 


0.062 


0.037 


0.043 


0.001 


0.002 


0.001 


Combinatoric bkg. 


0.019 


0.142 


0.045 


0.001 


0.005 


0.002 


Selection criteria 


0.068 


0.178 


0.046 


0.083 


0.172 


0.011 


Detector model 


0.064 


0.080 


0.064 


0.054 


0.040 


0.054 


Quadrature sum 


0.172 


0.251 


0.132 


0.122 


0.318 


0.083 



fects the level of mis-reconstructed events. 

To evaluate the effect of possible misalignments in the S VT 
on the mixing parameters, signal MC events are reconstructed 
with different alignment parameters, and the analysis is re- 
peated. The misalignments introduce resolution offsets in 
the MC of up to 10 fs and the corresponding fitted lifetimes 
change by up to 3 fs. However, since the lifetimes of all decay 
modes change by similar amounts, the effect on yep and AY 
is small. We also changed the energy loss correction applied 
in the tracking by 20% since a previous analysis has shown 
that the energy loss is underestimated in the reconstruction of 
data events 11311 . This changes the fitted lifetimes by about 
0.5 fs but has little effect on the mixing parameters. 

We combine the results shown in table [H with those from a 
previous BaBar study J3, based on 91 fb _1 of data, that 
does not require a D* + parent to identify the D° decays. 
While use of these untagged D° decays increases the sen- 
sitivity to yep through a factor of five increase in statistics, 
it also introduces different background behavior and there- 
fore different systematic errors. We have not used these un- 
tagged events in the current analysis, and thus the untagged 
data sample of the earlier analysis is essentially disjoint and 
its results statistically independent. Systematic uncertainties 
in the previous analysis were dominated by the limited num- 
ber of simulated events. Since the MC samples in the present 
study are presented here are entirely independent, this uncer- 
tainty is not correlated with those on the new results. Con- 
servatively assuming the remaining systematic uncertainties 
to be 100% correlated, we combine the two results using the 
BLUE method and obtain y C p = [1.03 ± 0.33(stat) ± 
0.19(syst))]%. 

In summary, we have obtained a value of yep = [1.24 ± 
0.39(stat) ± 0.13(syst))]% which is evidence of D°-D° 
mixing at the 3er level. It is compatible with our previous 
result lfl4ll and the recent lifetime ratio measurement from 
Belle of y C p = [1.31 ± 0.32(stat) ± 0.25(syst)]% 0]. We 
find no evidence for CP violation and determine AY to be 
[-0.26 ± 0.36(stat) ± 0.08(syst)]%. The result is consistent 
with SM estimates for mixing. 
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